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Abstract 

We are witnessing an era of intense experimental efforts that will provide in¬ 
formation about the properties of nuclei far from the valley of j5 stability, re¬ 
garding resonant and scattering states as well as (weakly) bound states. This 
talk describes our formalism for including these necessary ingredients into the 
No-Core Shell Model by using the Gamow Shell Model approach. Applica¬ 
tions of this new approach, known as the No-Core Gamow Shell Model, both 
to benchmai'k cases as well as to unstable nuclei will be given. 


1 Introduction 


The atomic nucleus, which at energies in the range of mega-electron volts can be viewed as a quantum 
system of strongly interacting protons and neutrons ffl] , is a very fascinating object. One of the most 
intriguing challenges is the description of its time evolution. Indeed the nucleus can be naturally found 
in quantum states that decay by emitting photons, electrons or even nucleons and heavier particles. This 
endeavor becomes even more necessary in the rise of new experiments that probe new modes of radioac¬ 
tivity which will need theoretical justification |0]. One way is to solve the time-dependent Schrddinger 
equation by means of time discretization techniques iM. A different approach lies in the solution of 
the time-independent Schrddinger equation in the complex energy or momentum space. The time depen¬ 
dence is then absorbed by the complex nature of the solution, whose imaginary part is associated with 
the decay time. Metastable nuclear states and resonances could then be described in a time-independent 
formalism imi. Formulating the structure and reaction problem in the complex energy plane 13- 
22n . provides with an alternative step towards the unification of structure (bound states) and reaction 
(resonances) aspects in nuclear physics, which will lead in a more controlled and model independent 
evaluation of observables. Noted that there already exists vivid progress on the ab-initio description of 
structure and reactions on the real-axis by Fawrence Fivermore and TRIUMF groups i23l425n. by Fos 
Alamos/Argonne groups |I3] and also on the lattice by HS- 


2 The NC(GSM) formalism 

One of the ways to obtain complex energy solutions of a physical system is by turning the Schrddinger 
equation into an eigenvalue problem and diagonalizing the complex Hamiltonian matrix. The basic code 
that has been employed for the description of resonances, by diagonalizing a very lyge non-Hermitian 
complex symmetric Hamiltonian matrix, is the Gamow Shell Model (GSM) code 129143311 . The orthonor¬ 
mal underlying basis, upon which the Hamiltonian matrix is diagonalized, is known as Berggren basis 
10] which provides a symmetric description of bound states, resonances and scattering states. The Hamil¬ 
tonian matrix is non-Hermitian and complex symmetric with complex eigenvalues. The matrix is sparse 
and a relatively small number of eigenvectors and eigenvalues is of interest. At this point there is no 
conceptual difference between GSM with a core and the NCGSM besides the underlying Hamiltonian. 
Indeed, in the GSM case the Hamiltonian consists of the one-body kinetic energy, the one-body mean- 
field (either schematic or Hartree-Fock) and the residual NN interaction, whereas in the NCGSM case the 
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one-body mean field is not present in the A-body Hamiltonian ll32ll . The nucleon-nucleon interaction is 
expressed in the Berggren basis using the potential separable expansion (PSE) 041] in a HO basis method 
1 32|, 133,13^- In this way, matrix elements between scattering and/or resonant states never diverge, due 
to the Gaussian fall-off of the HO radial form factor, and in addition one can conveniently transform 
matrix elements from the relative frame to the lab frame. For matrix elements of other operators, such 
as electromagnetic transition operators, the renormalization of integrals relies on the method of external 
complex scaling (ECS). We note here though that the PSE method has also been employed for the cal¬ 
culation of the recoil matrix elements of the intrinsic Hamiltonian, since in this case the ECS technique 
does not provide converged results. 

For the diagonalization of the large complex symmetric matrix we have used a complex extension 
of the Eanczos algorithm Bin . The largest matrix that has been diagonalized has a dimension of 10®. 
This number materializes to about A=4,5 for ab-initio no-core calculations in the Berggren basis and to 
7-8 active valence particles when assuming a frozen configuration (also known as “core") but allowing 
only a portion of the particles occupying the continuum (Berggren) orbitals (particle-hole truncations). 
The second alternative is the Davidson method Bill (see also a recent application on the description 
of tunneling for a two-body atomic system BTI] '). We would like to highlight a unique feature of the 
GSM diagonalization. Both the Eanczos and the Davidson methods calculate the ground state of the 
system as the lowest eigenvalue. In the GSM where the Hamiltonian matrix is non-Hermitian, the lowest 
eigenvalue is not guaranteed to coincide with the ground state and it may as well be a scattering state. A 
criterion is established which separates the state of interest from the wealth of scattering states; that is the 
overlap method 13111 . In the overlap method a smaller diagonalization takes place first, in a space spanned 
by a few states (usually single particle bound states and resonances). This smaller space sometimes is 


Table 1: Comparisons Il38l] between the NCGSM and the computer program MFDn. Ngheii = 2nH-f and hu = 20 
MeV. 


Nucleus 

MFDn 

NCGSM 

Difference 

^Hl+(N,hell = 4) 

-1.6284 

-1.6284 

< 0.1 keV 

2H 1+ (N.hell = 8) 

-2.1419 

-2.1419 

< 0.1 keV 

3H 1/2+ (N,hell = 4) 

-7.6016 

-7.6016 

< 0.1 keV 

3h 1/2+ (Nshell = 8) 

-8.3203 

-8.3203 

< 0.1 keV 

3He l/2+(N,heii = 8) 

-7.6084 

-7.6084 

< 0.1 keV 

%e 0+ (N.heii = 4) 

-27.3685 

-27.3684 

0.1 keV 

®Eil+(N,hell = 4) 

-24.9778 

-24.9776 

0.2 keV 

®Ei 3+ (Nsheii = 4) 

-22.4959 

-22.4957 

0.2 keV 


also called the “pole approximation". The solution is the reference eigenvector. At a second step a 
diagonalization in the full space (bound states-i-resonances-i-scattering states) takes place and the solution 
is the one that maximizes the overlap with the reference eigenvector. 


Together with the GSM, another algorithm known as the Density Matrix Renormalization Group 
(DMRG) has also been used 132 11391 - 14211 . The method is a truncation algorithm which aims in reducing 
the sizes of the matrices to be diagonalized, while keeping the same accuracy as in the full calculation. 
DMRG is an iterative method for which at each step the space is increased by adding basis states one- 
by-one and the truncation is dictated by the density matrix which is constructed in each step. At each 
step a complex symmetric non-Hermitian matrix is diagonalized via the Eanczos method and the matrix 
is smaller than the typical matrix of a full-fledged GSM calculation. Several diagonalizations of such 
smaller matrices are performed until a convergence criterion is reached. A recent variant of the DMRG 
method i43n has shown that very large model spaces can be reached and many “sweeps" can be per¬ 
formed in a timely manner, making it tempting to use such a variant also in Berggren basis DMRG based 
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calculations. 


3 Applications of the NCGSM approach 
3.1 Benchmarks 


Aiming at a predictive theory it becomes increasingly important to complete a quality control on the 
solvers which are employed for the solution of the many-body problem, a task that implies validation 
and cross-checks (benchmai'ks) of existing codes 14411 . To test the NCGSM algorithm we have performed 
benchmark calculations, in which we compared the NCGSM results with results obtained using the 
NCSM 1^1 Many-Fermion Dynamics nuclear (MFDn) computer program |4^, T7], and in Table ([T]) we 
present applications for systems up to ®Li. For this benchmark we employed a Harmonic Oscillator 
(HO) basis in a Full Configuration Interaction (FCI) truncation and the JISP16 realistic interaction 14811 . 
Nevertheless, chiral interactions were also tested successfully 0811 . This agreement reassures that the 
calculations are not contaminated with unintentional errors or flaws. It should be noted that there has 
also been two other successful works that benchmarked the GSM algorithm, using central interactions 
and a a-core, against the complex scaling technique |^, 


3.2 ANCs and widths 

In its current implementation the NCGSM is not ready to provide scattering observables on the real-axis, 
such as cross-sections, even though such a goal is not far from reach after the combination of the GSM 
with the Resonating Group Method o ill using phenomenological interactions/optical potentials 1151 15211. 
We are able however to compute overlaps between nuclear states and access information associated with 



Figure 1: (Color online) Overlap function and tail fit with a Hankel function. Cij stands for the ANC. Figure is 
fromlH. 


the “tail" of the overlap. The relevant quantity in the study of asymptotic properties of the nuclear wave- 
function or actually its projection onto cluster (sub-cluster) states 1153n is the Asymptotic Normalization 
Coefficient (ANC). Recently there is a collective effort in nuclear theory to compute asymptotic quanti¬ 
ties and we are witnessing an abandonment of quantities such as spectroscopic factors in favor of ANCs 
and widths or partial widths. The basic argument behind this endeavor, besides the physics interest (e.g. 
relevance to astrophysical processes for both resonance widths and ANCs), is the fact that asymptotic 
quantities are less model dependent and closer to the notion of an observable quantity 1541. 15511. At the 
same time ANCs can serve as an internal consistency test between many-body methods since calculations 
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at distances far away from the nuclear interaction range always pose challenges and difficulties (see e.g. 
discussion at 11561] for some of the ab-initio methods). The GSM or the NCGSM which are formulated 
on a basis that has a correct asymptotic behavior and captures the relevant long range physics, become 
appropriate for the calculation of asymptotic quantities. For a detailed review of the progress that has 
been made in the calculation of ANCs and also the experimental situation we refer the reader to 115611 and 
also 157] for ANCs calculations within the GSM. 


In this contribution we present calculations that were published in ll32|] of ANCs within the 
NCGSM using realistic interactions. The model space of the calculation, as it was described in i32ll . 
includes single particle (s.p.) partial waves with angular momentum up to ^ = 4 (g-waves). For ®He, 
being particle unstable in its ground state, we employed a complex GHF basis consisting of the s.p. 
0 p 3/2 resonant state and non-resonant P 3/2 states along the complex contour which encloses the s.p. res¬ 
onant state, a necessary condition for the s.p. Berggren completeness to be satisfied. In order to obtain 
the many-body solution and calculate the overlap for the reaction: ^Heo 4 . + ^^ 63 / 2 - we used the 

Davidson diagonalization 13IH method and limited our selves to continuum configurations that allowed 
up to four particles occupying continuum orbits (4p-4h truncation). In the future it will be important to 
assess the importance of the configurations involving many particles in continuum orbits and also try to 
accompany the results with a truncation error associated with the missing configurations. It should also 
be noted that the creation operator for the calculation of the overlap does not only create on a single state, 
but there is a sum over all continuum states; a fundamental difference between NC(GSM) and traditional 
configuration interaction calculations in a HO basis Is^ . 

We present our results in Fig U] wh ich was taken from |32] for the effective yiowk A = 

1.9fm~^ chiral N^LO NN interaction 16211 . We see that the overlap exhibits both real and imaginary parts 
reflecting the complex nature of the Berggren basis. After fitting the asymptotic part of the overlap with 
a complex Hankel function we extract an ANC with a real part of 0.197 fm“^/^. Now knowing the ANC 
we are able to obtain the width of ^He using the formula: 


C = 


Tfi 


hm{k) 


( 1 ) 


that relates the width and the ANC 115TP . where // is the reduced mass, k is the real part of the complex 
linear momentum that corresponds to the neutron-separation energy of ®He and C stands for the ANC. 
The result for the width is T = 311 keV. Within the 4p-4h truncation, the complex Davidson diagonal¬ 
ization provided for the Si,i the value of -1.561 MeV and a width T of 370 keV. The small difference 
between the width obtained from the ANC formula and the one obtained from the diagonalization stems 
from an approximation that was made on the proof of formula ([Til. The approximation is that the real 
part of the linear momentum has to be considered. It has been shown in i32ll that this approximation 

implies the condition - 0 , namely, formula ([Hi will work for states that have widths much smaller 

as compared to two-times the separation energy. In our case = 11.8 %, which explains the deviation 
from the exact diagonalization result. 


3.3 NCGSM results for and ^Li 

Next we are presenting calculations for and ^Li. These systems have also been computed in ap¬ 


proaches with realistic interactions, in Il63[] for a scattering benchmark calculation within several few 
body methods and in 16411 within the complex scaling method. Especially for hydrogen isotopes there 
exists a recent experimental interest since measurements claim to observe a relatively narrow res¬ 
onance ll65ll . In addition, it is believed that the triplet of isotopes resemble the Helium isotopic 

chain (®’^’®He). Namely ^H resembling ®He is the most “bound" member of the triplet having a relatively 
small width, ®H resembling ^He is “unbound" reflecting its extremely large width while ®H is just a broad 
resonance (broader than ^H), resembling in this case ®He which is less bound than ®He. The theoretical 
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Figure 2: (Color online) Ground state and first excited state spectra of and ^Li. Gray scale denotes the 
continuum/scattering regime which is described by adopting a complex Beggren basis. 


investigation of the shell structure and the pairing correlations in the continuum will be a decisive factor 
in the understanding of the binding mechanism in this area of the chart and it will also be a challenging 
task. 

In this work we are using a phenomenological WS potential for the generation of the basis. For '^H 
and ^Li the WS basis is created for the + n or ^He + p systems respectively. Namely for a neutron 
0 p 3/2 s.p. resonant state is considered whereas for ^Li a proton OP 3/2 s.p. resonant state. The long-range 
Coulomb interaction in the case of ^Li is treated in the same way as in 11321. rod l67n. Besides the complex 
neutron and proton resonances and the associated complex continua, we consider s.p. partial waves up to 
h-waves (^=5) as HO basis states characterized by a frequency of /icc = 20 MeV. In particular, for and 
for the neutron space we used a basis that consisted of 15p3/2 complex s.p. states and 1 ^ 1 / 2 , 7pi/2, 2d, 
2f, 2g, Ih real HO s.p. states, while for the proton space all basis states are HO states with 7si/2! ^Vi/ 2 ^ 
2d, 2f, 2g, Ih. Noted that we are using a different truncation scheme that departs from the typical N^ax 
or Ngheii HO truncation. Here, and also in other NCGSM calculations, we choose to use more radial 
nodes for specific £ states. Here states with £=0 ,1 have been chosen to have more nodes since we noticed 
that they contribute more energy to the system. Further calculations that investigate this behavior are in 
progress. The same basis was used for ^Li but instead of a neutron s.p. resonant state we considered 
a proton s.p. resonant state. Our results for the g.s and first excited states of ^H and “^Li are shown 
on FigjJl The energies are with respect to the one-neutron and one-proton particle thresholds. For the 
thresholds, the g.s. energies of ^H and ^He are found to be -7.92 MeV and -7.12 MeV respectively, for 
an effective Viowk ^ = 2.0fm“^ N^LO Chiral EFT interaction. Overall we observe a good agreement 
with experimental measurements llhSH and especially the small gap between the 1” and 2~ states that 
has been observed is also predicted in our calculation. Even though the calculations for the widths show 
stability, the widths are very large, so we are not dealing with typical resonances that would have an 
appreciable impact on cross-sections. In the calculation presented here we restricted the occupation of 
continuum orbits to a maximum of three particles (3p-3h). In our future work we will provide results 
without truncations and also with investigations on the impact of missing excitations on energies and 
widths. 


4 Conclusions and outlook 

In conclusion, we presented applications of the NCGSM for the calculation of energies, widths and 
asymptotic quantities such as ANCs of unbound nuclei. At the same time we benchmarked our algorithm 
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against another commonly used solver. Our immediate goals are to continue working on applications 
of the NCGSM for light unbound nuclear systems in the neutron and proton rich side of the nuclear 
chart. The success of the NCGSM is tied with advances in computer algorithms and the formulation of 
efficient complex symmetric diagonalization solvers. On the physics side, even when the most efficient 
solver will be at hand, a full calculation in the Berggren basis (no p-h truncations) for a system such 
as ^^Li will be very demanding. Hence, one of our goals is to quantify at first the impact of missing 
truncations (we already know that the weight of many particles in continuum orbits is small). Also 
we aim at a construction of a realistic effective interaction, i.e. an interaction in the NCGSM nuclear 
medium Q that would stem from a realistic free-space interaction, which we will then utilize for NCGSM 
calculations (see e.g. Ihil] for a Lee-Suzuki transformation for interactions in the complex energy plane 
and a multireference perturbation theory approach). Finally, we would like to contribute to the effort 
of bridging the gap between Lattice-QCD and low energy physics and a possible way was shown in 
117on . This could be achieved once NCGSM will be utilized to handle NNN interactions. Then we could 
naturally compute resonant features of finite nuclei at any pion mass that is available at that time. 
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* NCGSM medium here should be seen as the model space that the A-nucleons evolve, but in addition, scattering to contin¬ 
uum orbits will be allowed through the Berggren basis 
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